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Abstract 

We analyse some peculiar properties of the function of the Mittag-Leffler 
(M-L) type, ea{t) := Ea{—t") for < a < 1 and t > 0, which is known to 
be completely monotone (CM) with a non negative spectrum of frequencies 
and times, suitable to model fractional relaxation processes. We first note 
that (surprisingly) these two spectra coincide so providing a universal scaling 
property of this function, not well pointed out in the literature. Furthermore, 
we consider the problem of approximating our M-L function with simpler CM 
functions for small and large times. We provide two diff'erent sets of elementary 
CM functions that are asymptotically equivalent to ea{t) as t — )• and t — t- -|-oo. 
The first set is given by the stretched exponential for small times and the 
power law for large times following a standard approach. For the second set 
we chose two rational CM functions in t", obtained as the Fade Approximants 
(FA) [0/1] to the convergent series in positive powers (as t — t- 0) and to the 
asymptotic series in negative powers (as t — t- do), respectively. From numerical 
computations we are allowed to the conjecture that the second set provides 
upper and lower bounds to the Mittag-Leffler function. 

1 Introduction 

Since a few decades the special transcendental function known as Mittag-Leffler func- 
tion has attracted an increasing attention of researchers because of its key role in 
treating problems related to integral and differential equations of fractional order. 



Since its introduction by the Swedish mathematician Mittag-Leffler at the beginning 
of the last century up to the 1990's, this function was seldom considered by mathe- 
maticians and applied scientists. 

Prom a mathematical point of view we recall the paper by Hille and Tamarkin (1930) 
on the solutions of the Abel-Integral equation of the second kind, the Handbook of 

the Bateman project, Erdelyi (1955), where it was restricted in the chapter devoted 
to miscellaneous functions, and the books by Davis (1936), by Sansone and Gcrretsen 
(1960), by Dzherbashyan (1966) (unforunately in Russian) and finally by Samko, 
Kilbas and Marichev (1987-1993). 

For apphcations we recall a few papers by Cole (1933), Cole and Cole (1941, 1942), 
Gross (1947) and Caputo and Mainardi (1971a,1971b) where the Mittag-Leffler func- 
tion was adopted to represent the responses in dielectric and viscoelastic media. 

In the 1960's the Mittag-Leffler function started to exit from the realm of miscella- 
neous functions because it was considered as a special case of the general class of Fox 

H functions, that can exhibit an arbitrary number of parameters in their integral 
Mellin-Barnes representation, see e.g. the books by Kiryakova (1994), Kilbas and 
Saigo (2004), Mathai et al. (2010). 

However, in our opinion, this classification in a more general framework has, to some 
extent, obscured the relevance and the applicability of this function in applied science. 
In fact most mathematical models are based on a small number of parameters, say 
1 or 2 or 3, so that a general theory may be confusing whereas the adoption of 
a generalized Mittag-Leffler function with 2 or 3 indices may be sufficient, see e.g. 
Beghin and Orsingher (2010), Tomovski et al. (2010), Capelas et al. (2011), Sandev 
et al. (2012). Multi-index Mittag-Leffier functions have been introduced as well, see 
e.g. Kiryakova (2010), Kiryakova and Luchko (2010), Kilbas et al. (2013), but their 
extensive use has not yet been pointed out in applied sciences. 

Nowadays it is well recognized that the Mittag-Leffler function plays a fundamental 
role in Fractional Calculus even if with a single parameter (as originally introduced 
by Mittag-Leffler) just to be worth of being referred to as the Queen Function of 
Fractional Calculus^ see Mainardi and Gorcnfio (2007). On this respect we just point 
out a few of reviews and treatises on Fractional Calculus (in order of publication 
time): Gorenflo and Mainardi (1997), Podlubny (1999), Hilfer (2000), Kilbas et al. 
(2006), Magin (2006), Mathai and Haubold (2008), Mainardi (2010), Diethelm (2010), 
Klafter et al (2012), Baleanu et al (2012). 
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2 The Mittag-Lefler function in fractional 
relaxation processes 

The Mittag-Leffler function is defined by the following power series, convergent in the 
whole complex plane, 

oo n 

Wc recognize that it is an entire function providing a simple generalization of the 
exponential function to which it reduces for a = 1. We also note that for the con- 
vergence of the power series in (2.1) the parameter a may be complex provided that 
3?(q;) > 0. The most interesting properties of the Mittag-Leffler function are asso- 
ciated with its asymptotic expansions as 2; — )■ oo in various sectors of the complex 
plane. For detailed asymptotic analysis, which includes the smooth transition across 
the Stokes lines, the interested reader is referred to Wong and Zhao (2002). 

In this paper we limit ourselves to the Mittag-Leffler function of order a G (0,1) on 
the negative real semi- axis where is known to be completely monotone (CM) due a 
classical result by Pollard (1948), who proved a conjecture by Feller. 

Let us recall that a function (j){t) with t G is called a completely monotone (CM) 
function if it is non negative, of class C°°, and {—l)^(j)^"-\t) > for all n G IN. Then 
a function ip{t) with t G lF?t is called a Bernstein function if it is non negative, of 
class C°°, with a CM first derivative. These functions play fundamental roles in linear 
hereditary mechanics to represent relaxation and creep processes, see e.g Mainardi 
(2010). For mathematical details we refer the interested reader to the survey paper 
by Miller and Samko (2001) and to the recent book by SchiUing et al. (2010). 

In particular we are interested to the function 

00 j.an 

e^{t):^E^{-r) = Y.{-ir——- t>0, 0<a<l, (2.2) 

n=0 ^ l*^^ ^) 

that provides the solution to the fractional relaxation equation, see Gorenflo and 
Mainardi (1997), Mainardi and Gorenflo (2007), Mainardi (2010). 

For readers' convenience let us briefly outline the topic concerning the generaliza- 
tion via fractional calculus of the first-order differential equation governing the phe- 
nomenon of (exponential) relaxation. Recalling (in non-dimensional units) the initial 
value problem 

djU 

— = -u(t) , t > , with m(0+) = 1 (2.3) 
at 
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whose solution is 

u{t) — exp(— t) , 

the following two alternatives with a e (0, 1) are offered in the literature: 



du 
'di 



D]-"" u{t) , t>0, with u{0+) = 1 , 



u{t) = -u{t) , t>0, with u{0' 



(2.4) 



(2.5a) 
(2.56) 



where -D^ ~" and *Z^" denote the fractional derivative of order 1 — « in the Riemann- 
Liouville sense and the fractional derivative of order a in the Caputo sense, respec- 
tively. 

For a generic order /i e (0, 1) and for a sufficiently well-behaved function f{t) with t G 
IRt the above derivatives are defined as follows, see e.g. Gorenflo and Mainardi(1997), 
Podlubny (1999), 



1 



d 



r(i - //) dt 
1 



fir) 



{t- tY 

fir 



dr 



r(l - /i) Jo {t - tY 
Between the two derivatives we have the relationship 



r(i-/.) 



dr . 



fit) - fio' 



(2.6a) 
(2.66) 

(2.7) 



Both derivatives in the limit // — > 1 reduce to the standard first derivative but for 
— >■ 0+ we have 



D!;f{t)^f{t), Mfit) = fit)-fio^), /^^o+. 



(2.^ 



In analogy to the standard problem (2.3), we solve the problems (2.5a) and (2.5b) 
with the Laplace transform technique, using the rules pertinent to the corresponding 
fractional derivatives. The problems (a) and (b) are equivalent since the Laplace 
transform of the solution in both cases comes out as 

-(^) - ^ ' (2-9) 

that yields our function 

u{t) = eait) ■.= E^{-e). (2.10) 

The Laplace transform pair 

e„(t)-^^, a>0, (2.11) 

S ~r J- 



4 



can be proved by transforming term by term the power series representation of ea{t) 
in the R.H.S of (2.2). Furthermore, by anti-transforming the R.H.S of (2.10) by 
using the complex Bromwich formula, and taking into account for < o; < 1 the 
contribution from branch cut on the negative real semi-axis, we get 



POO 

Jo 

where, for the Titchmarsh formula. 







1 


^ - Im < 






TT 




s=r e^*"^ y 



r" ^ sin (an) 



+ 2 r"' cos (avr) + 1 



> 0. 



(2.12) 



(2.13) 



Since K^ir) is non negative for all r in the integral, the above formula proves that 
ea{t) is CM function in view of the Bernstein theorem. This theorem provides a 
necessary and sufficient condition for a CM function as a real Laplace transform of a 
non negative measure. 

However, the CM property of ea{t) can also be seen as a consequence of the result 
by Pollard because the transformation x — t"' is a, Bernstein function for a e (0, 1). 
In fact it is known that a CM function can be obtained by composing a CM with a 
Bcrnstrcin function based on the following theorem: Let (f){t) be a CM function and 
let ip{t) be a Bernstein function, then 4>[ip{t)] is a CM function. 

As a matter of fact, Ka{r) provides an interesting spectral representation of ea{t) in 
frequencies. With the change of variable t — 1/r we get the corresponding spectral 
representation in relaxation times, namely 



JO 



(2.14) 



that can be interpreted as a continuous distributions of elementary (i.e. exponential) 
relaxation processes. We also note that (surprisingly) we get the identity between the 
two spectral distributions 

K^{r) = H^{T), (2.15) 

a fact pointed out in Linear Viscoelasticity by the author in his recent book, Mainardi 
(2010). This kind of universal/scaling property seems a peculiar one for our Mittag- 
Leffier function ea{t). In Fig 1 we show Ka{r) for some values of the parameter a. 
Of course for a = 1 the Mittag-Leffler function reduces to the exponential function 
exp(— t) and the corresponding spectral distribution is the Dirac delta generalized 
function centred at r = 1, namely 6{r — 1). 
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Fig.l Plots of the spectral function Kaijr) for a = 0.25,0.50,0.75,0.90 in the fre- 
quency range < r < 2. 



3 Asymptotic approximations to the Mittag-Lefler 
function 



In Fig 2 we show some plots of ea{t) for some values of the parameter a. It is worth 
to note the different rates of decay of ea{t) for small and large times. In fact the 
decay is very fast as t — )■ 0+ and very slow as t — +oo. 




Fig.2 Plots of the Mittag-Leffler function e„(t) for a = 0.25, 0.50, 0.75, 0.90, 1. in the 
time range < t < 15 



3.1 The two common asymptotic approximations 

It is common to point out that the function ea{t) matches for t — t- 0^ with a stretched 
exponential with an infinite negative derivative, whereas as t — )■ oo with a negative 
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power law. The short time approximation is derived from the convergent power series 
representation (2.2). In fact, 



ea{t) = 1 



r(l + a) 



+ . . . ~ exp 



r(l + a) 



(3.1) 



The long time approximation is derived from the asymptotic power series represen- 
tation of ea{t) that turns out to be, see Erdelyi (1955), 



ea(i)-^E(-l) 



n-1 



so that, at the first order. 



n=l 



ea{t) 



r(l — an) 



i — >■ oo , 



r(i - a) ' 



t ^ oo . 



(3.2) 



(3.3) 



It is common to say that the function ea{t) interpolates the stretched exponential due 
its very fast decay for short times, whereas interpolates the negative power law due to 
its very slow decay for long times. Then, wc get the two asymptotic representations 
as t ^ and t ^ oo commonly found in the literature: 



ea{t) ^ < 



e° (t) := exp 



t~°' sin(Q;7r) T{a) 



ni-a) 



TT 



i — >■ oo . 



(3.4) 



The stretched exponential has been adopted in order to replace the divergent expres- 
sion 1 — t"/r(l -|- a) obtained from (3.1), asymptotically equivalent to ea{t) as i — >■ 0. 
Of course we have the inequality 



elit)<e^{t), t>0, 0<a<l. 



(3.5) 



In Figs 3-7 LEFT we compare for a = 0.25, 0.5, 0.75, 0.90, 0.99 in logarithmic scales 
the function ea{t) and its asymptotic representations, the stretched exponential ej^(t) 
valid for t — ?> (dashed line) and the power law e^(t) valid for t — ?> cxo (dotted line). 
We have chosen the time range 10^^ < ^ < lO"*"^. In the RIGHT we have shown the 
plots of the relative errors (in absolute values) 

K{t)-e^{t)\ |e-(t)-e.(t)| 

pointing out a line at an error 1% under which the approximations can be considered 
reliable. 
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Fig. 3 Approximations e°(t) (dashed line) and e^(t) (dotted line) to Cait) (LEFT) 
and the corresponding relative errors (RIGHT) in 10~^ < i < lO"*"^ for a = 0.25. 




Fig. 4 Approximations e°(t) (dashed line) and e^(t) (dotted line) to ea{t) (LEFT) 
and the corresponding relative errors (RIGHT) in 10^^ < t < 10^^ for a = 0.50. 




Fig. 5 Approximations e°(t) (dashed line) and e^(t) (dotted line) to ea(t) (LEFT) 
and the corresponding relative errors (RIGHT) in 10~^ < t < 10^^ for a = 0.75. 
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Fig. 6 Approximations e°(t) (dashed line) and e^(t) (dotted line) to Cait) (LEFT) 
and the corresponding relative errors (RIGHT) in 10^^ < t < 10^^ for a = 0.90. 




Fig. 7 Approximations e°(t) (dashed line) and e^(t) (dotted line) to ea(t) (LEFT) 
and the corresponding relative errors (RIGHT) in 10^^ < ^ < 10^^ for a = 0.99. 

We note from Figs 3-7 that, whereas the plots of e^{t) remain always under the 
corresponding ones of ea{t), the plots of e^(t) start upper those of Cait) but, at a 
certain point, an intersection may occur so changing the sign of the relative errors. 



3.2 The two rational asymptotic approximations 

We now propose a new set of CM functions approximating Cait): {fa(t), Qait)}, alter- 
native to {e° (t), e^(^)}, obtained as the first Fade approximants [0/1] to the power 
series in (2.2) and (3.2), respectively. For more details on the theory of Fade 
Approximants we refer e.g. to Baker (1975). We thus obtain the following rational 
functions in 

fait):=—^^r^l-—^—^e^it), t^O, (3.7) 

i + f(iT^o r(i + «) 
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9ait) := , ^.},^ r ~ ^..^ " . ~ e„(t) ,t^oo. (3.8) 

l+t°r(l — a) r(i — a) 

Now we prove the inequality 

ga{t)<Ut), t>0, 0<a<l, (3.9) 

as a straightforward consequence of the reflection formula of the gamma function. In 
fact, recalling the definitions (3.7)-(3.8) we have for Vt > and a e (0, 1): 

1 Tia 
9a{t)<Ut) ^ r(i-a)>— — ^ r(l-a)r(l + a) = — ^>1. 

1 ( 1 + a) smi vra) 



In Figs 8-12 LEFT we compare for a = 0.25, 0.5, 0.75, 0.90, 0.99 in logarithmic scales 
the function ea{t) and its rational asymptotic representations, fa{t) valid for t — )■ 
(dashed line) and ga{t) valid for t — )■ oo (dotted line). We have chosen the time range 
10~^ < t < lO"*"^. In the RIGHT we have shown the plots of the relative errors (no 
longer in absolute values) 

fa{t)-ea{t) ea{t)-g^{t\ 

pointing out a line at an error 1% under which the approximation can be considered 
reliable. 




Fig. 8 Approximations /^(t) (dashed line) and Qaif) (dotted line) to Caif) (LEFT) 
and the corresponding relative errors (RIGHT) in 10~^ < t < lO"*"^ for a = 0.25. 
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Fig. 9 Approximations fa(t) (dashed line) and ga(t) (dotted line) to Cait) (LEFT) 
and the corresponding relative errors (RIGHT) in 10~^ < ^ < lO"*"^ for a = 0.50. 




Fig. 10 Approximations fait) (dashed line) and gait) (dotted line) to Cait) (LEFT) 
and the corresponding relative errors (RIGHT) in 10~^ < t < lO"*"^ for a = 0.75. 




Fig. 11 Approximations fait) (dashed line) and gait) (dotted line) to Cait) (LEFT) 
and the corresponding relative errors (RIGHT) in 10^^ < t < 10^^ for a = 0.90. 
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— e„(t) = E.(-f°) 

--/.(*)= [i + nfc]]"' 
■■s„(t) = [i + t»r(i-a)]-' 




--!/„(() -£'„Wl/e.(f): /o(t)= [l + r(^J 
■■ [e„(t) -9„(f)l/e.W: 94f) = |1 +ff(l -q)!" 



a = 0.99 



Fig. 12 Approximations fait) (dashed line) and gait) (dotted line) to Cait) (LEFT) 
and the corresponding relative errors (RIGHT) in 10~^ < t < lO"*"^ for a = 0.99. 



As a matter of fact, from the plots in Figs 8-12, we recognize, for the time range and 
for a E (0, 1) considered by us, the relevant inequality 

9ait) < eait) < fait) , (3.11) 

that is gait) and fait) provide lower and upper bounds to our Mittag-LefHer function 
Cait) . This of course can be seen as a conjecture that we leave as an open problem 
to be proved (or disproved) by specialists of CM functions. 

We also see that whereas the short time approximation fait) turns out to be good 
only for small times, the long time approximation Qait) is good (surprisingly) also 
for short times, falling down only in an intermediate time range. In fact from the 
RIGHT of Figs 8-12 we can estimate the ranges of validity when the relative error is 
less than 1%. We could further discuss on them. 

Concerning Fade Approximants (PA), we like to cite the paper by Freed et al. (2002), 
where an appendix is devoted to the table of PA for the Mittag-Leffler function 
Eai—x), a; > 0. We point out the fact that our first PA's [0/1] are constructed from 
two different series (in positive and negative powers of t") both resulting CM functions 
whereas the successive PA's of higher order can exhibit oscillations for t > 0. 



4 Conclusions 

We have discussed some noteworthy properties of the Mittag-Leffler function Eai—t"") 
with < a < 1 in the range t > 0. 

Being a completely monotone (CM) function, because of the Bermstein theorem this 
function admits non negative frequency and time spectra. We have pointed out that 
these two spectra are equal, so providing a universal scaling property. 
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Furthermore, in view of a numerical approximation we have compared two different 
sets of approximating CM functions, asymptotically equivalent to Ea{—t°') for small 
and large times: the former commonly used in the literature, the latter probably 
new. The last set is constituted by two simple rational functions that provide upper 
and lower bounds, at least in our numerical examples on a large time range. We are 
allowed to the conjecture that this bounding property is always valid for any t > 
and for any a G (0, 1): an open problem offered to specialists of CM functions. 

Last but not the least, we note that our rational function approximating Ea{—t'^) 
for large times provides (surprisingly) a good numerical approximation also for short 
times, so failing only in an intermediate range. 
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